Giant enhancement of photodissociation of polar dimers in electric fields 
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We explore the photodissociation of polar dimers in static electric fields in the cold regime using 
the example of the LiCs molecule. A giant enhancement of the differential cross section is found for 
laboratory electric field strengths, and analyzed with varying rovibrational bound states, continuum 
energies as well as field strengths. 



I. INTRODUCTION 

The availability of cold to ultracold polar dimers in 
their absolute ground state [1-3] represents a break- 
through towards the control of all molecular degrees of 
freedom, i.e., the center of mass, electronic, rotational 
and vibrational motions. Such a polar sample provides 
a unique starting point to manipulate the internal struc- 
ture of molecules, and therefore the molecular dynamics 
with electrostatic fields. Oriented pendular states present 
a paradigm in this context. They exhibit a preferred spa- 
tial direction of their permanent electric dipole moments, 
which introduces an anisotropic and long-range electric 
dipole-dipole interaction between two dimers. Thus, by 
changing the directional features one can control and ma- 
nipulate physical processes taking place in the ultracold 
regime, e.g . , scattering dynamics 0-0 > chemical reac- 
tions spin dynamics [l2|,or even more, create 
quantum computational devices [13l4l5|. 

Here, we investigate a specific process, namely, the dis- 
sociation of polar diatomic molecules into two cold atoms 
via single-photon absorption, and how this process can 
be controlled and manipulated by using an additional 
external electric field. As a representative example, we 
consider the electronic ground state of the LiCs dimer 
and analyze the direct dissociation within the electronic 
ground state. The dependence of the angular distribu- 
tion of the photofragments on the initial rovibrational 
bound state, on the energy of the continuum, and on 
the static electric field strength are analyzed in detail. 
Our focus is the cold regime, where transitions to a few 
scattering partial waves contribute. In the later case the 
spectrum shows shape resonances showing up in the fact 
that, e.g., an initial bound state with rotational quan- 
tum number J ± 1 leads to large field-free differential 
cross sections and angular distributions via a shape res- 
onance of angular momentum J. Switching on a static 
field, the hybridization of the angular motion, i.e., mixing 
of the rotational states, takes places, and the dissociation 
process is significantly modified. We demonstrate that a 
moderate static electric field provokes enhancements of 
several orders of magnitude of the angular distributions 
of the dissociation process. 

The paper is organized as follows. In Sec. [Ill the rota- 
tional Hamiltonian is presented together with the differ- 



ential cross section and angular distribution of a direct 
dissociation process within the electronic ground state. 
We discuss the numerical results for LiCs in Sec. HIH 
analyzing the field- free and field-dressed cases. The con- 
clusions and outlook are provided in Sec. IIVI 

II. THE ROVIBRATIONAL HAMILTONIAN 
AND THE DISSOCIATION PROCESS 

The description of this process is achieved within the 
Born-Oppenheimer separation of the electronic and nu- 
clear coordinates. We assume that for the considered 
regime of field strengths a perturbation theoretical treat- 
ment holds for the description of the interaction of the 
field with the electrons, whereas a nonperturbative treat- 
ment is indispensable for the corresponding nuclear dy- 
namics. Thus, the rovibrational Hamiltonian of a polar 
dimer exposed to an homogeneous and static electric field 
reads 

H = T R + + V(R) - FD(R) cos 9 (1) 

where R and qb are the inter nuclear distance and the 
Euler angles, respectively, and we use the molecule fixed 
frame with the coordinate origin at the center of mass of 
the nuclei. Tr, hJ(9,(j)), and V(R) are the vibrational 
kinetic energy, orbital angular momentum, reduced mass 
of the nuclei and the field-free adiabatic electronic poten- 
tial energy curve, respectively. The last term in Eq. (pp) 
provides the interaction between the electric field and 
the molecule via its permanent electronic dipole moment 
function D(R). The electric field is oriented along the z- 
axis of the laboratory frame and posses the strength F. 
Our study is restricted to a non-relativistic treatment 
and addresses the spin singlet electronic ground state. 

In field-free space, each bound state of the molecule 
is characterized by its vibrational, rotational, and mag- 
netic quantum numbers (^, J, M). In the presence of the 
electric field only the magnetic quantum number M is 
conserved. However, for reasons of addressability we will 
refer to the electrically dressed states by means of the 
corresponding field- free quantum numbers (V, J, M). 

We assume a direct dissociation process within the 
electronic ground state with X E + symmetry caused by 
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a beam of linearly polarized light. The polarization of 
this laser is taken parallel to the direction of the static 
field. The initial state consists of a rovibrational bound 
state of the polar dimer exposed to an electrostatic elec- 
tric field and via the absorption of a photon a transition 
takes place to a continuum final state. The absorption 
of the photon from a linear polarized laser is described 
in terms of perturbation theory. Thus, making use of 
the dipole approximation the differential cross section is 
given by 

with = /(^; 77, E, F) being the angular distribution 

1(0) = \(9 c (n;E)\D(R) cos0|ig| 2 , (3) 

where ^ = (i?, 0, </>) is the rovibrational bound eigen- 
state of the Hamiltonian (pQ). It is characterized by the 
quantum numbers r] = (v, J, M), i.e., v and J are the 
field-free quantum numbers of this field-dressed level. 
The final state ^ c (tt;E) = ^ c (Q,R,0,^,E) is, in the 
field- free case, given by [l6j 

*c(0; E)=Y^ i 2J ' + l)e iA ""' i)j>M> (R, 0, 0; E) 

J'M> 

xD J M , o (n,0) (4) 

where i/jj'm'(R,6,(I>i E) is the energy normalized eigen- 
function of the Hamiltonian (pQ) describing the relative 
motion of the nuclei with kinetic energy E = hk 2 /2/z > 0. 
It is characterized by the phase 5j>m' = Aj>m' + 
through the relation ipj'M'(R ~^ 0? & E) ~ sm(kR + 
Sj'M' ~ kJ' /2)Yj/ M '(0,0), where Yj'M'(0,</>) are the 
spherical harmonics normalized-to-unity. The momen- 
tum k, with Euler angles Vt = (0,$), represents the 
propagating vector pointing along the final recoil direc- 
tion of the photofragments. D^, (Q, 0) = D^, (Q, <£>, 0) 
is a Wigner rotation matrix jl7j ]. 

We remark that if we assume the molecular ensem- 
ble to be randomly oriented, the population is equally 
distributed among the states with different values of M, 
and using the axial recoil approximation [l6| , the angular 
distribution (j3]) is reduced to the well known expression 
7(9, $) = cr(l+/3P 2 (©))/47r, with /3 being the anisotropy 
parameter, ^(0) the second-order Legendre polynomial, 
and it is normalized to the cross section a. This expres- 
sion was derived in the pioneering study of R.N. Zare and 
D. R. Herschbach [18], which was followed by a vast and 
rich amount of theoretical and experimental works. In 
particular, several of them have been devoted to the in- 
vestigation of the photodissociation of oriented molecules 

Since the probe laser is linearly polarized, we encounter 
in the absence of the static electric field the selection rules 
are A J = ±1 and AM = 0. Thus, the angular distribu- 
tion (j3j) is reduced to the square of a combination of two 
Legendre polynomials of degree J + 1 and J — 1, and it 



is symmetric with respect to 6 = tt/2. By turning on 
the electrostatic field, the hybridization of the rotational 
motion takes place and only the selection rule AM = 
holds. We restrict our analysis to initial and final states 
with M = 0, and due to the azimuthal symmetry the 
angular distribution is independent of the Euler angle 
<£, i.e., J(O) = 1(0). In our description, we consider 
the impact of the electric field on both the continuum 
and bound levels, and their wave functions are obtained 
from solving the Schrodinger equation associated to the 
rovibrational Hamiltonian (pQ). However, due to the long- 
range behaviour of the dipole moment, D(R) — >• D^R~ 7 , 
with D7 being constant, for R — >• 00, the field influence on 
the continuum states is smaller compared to the bound 
ones. Thus, we assume that there is no appreciable hy- 
bridization of its angular motion, so that the field-free 
rotational quantum number can be identified, and the 
expression for the wave function (j4]) is still holds approx- 
imately. 

Our aim is to explore the properties of the angular dis- 
tribution as the initial state, the continuum energy of the 
photofragments, and the field strength are varied. The 
two-dimensional Schrodinger equation associated with 
the nuclear Hamiltonian (pQ) has been numerically solved 
by employing a hybrid computational method that com- 
bines discrete and basis-set techniques applied to the ra- 
dial and angular coordinates, respectively [13, HH. The 
continuum is discretized and the numerically obtained 
wave functions are L 2 -normalized. However, the con- 
tinuum wave functions appearing in Eqs. (j3j) and ((4]) are 
energy-normalized. We address this issue with the help of 
a time-dependent formalism [27j. The continuum phase 
is computed by using the fact that our numerical tech- 
nique selects only those continuum states that are zero 
on the borders of the discretization box. 



III. RESULTS 

We investigate the impact of a static electric field on 
the dissociation of a polar dimer. As prototype example 
we take the 7 Li 133 Cs dimer in its electronic ground state 
X X E + . The potential energy curve is obtained from ex- 
periment [29|| including the correct long-range behavior. 
For the electric dipole moment function, we have used 
semiempirical data [30| , which were linearly extrapolated 
to short distances, and at long distance we have fitted 
them to the asymptotic behaviour D7/R 7 . Since D(R) 
is negative, in a strong electric field the pendular levels 
are antioriented along the field direction [27|, l3lM34j . The 
electronic ground state of 7 Li 133 Cs accommodates 55 vi- 
brational bands, and the last one v = 54 has two bound 
rotational levels with J = and 1 [32]. The contin- 
uum spectrum presents a shape resonance with energy 
E w 1.36 mK and rotational quantum number J = 2. 
Angular distribution will be provided in atomic units. 

Let us first analyze the dissociation process in the ab- 
sence of the electrostatic field. Figures QJa) and (b) 
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Fig. 1. Field- free angular distribution for the dissociation process 
to a continuum state with energy E ~ 0.3 mK, the initial states 
are (45, J, 0) (a) and (50, J, 0)(b) with rotational quantum numbers 
J = (full line),l (dashed line) and 2 (dotted line). 



present the angular distribution as a function of the angle 
6 for the initial states with v = 45 and 50 and J = 0, 1, 
and 2, respectively, and the final continuum energy is 
0.3 mK. For a certain rotational quantum number, /(©) 
is symmetric with respect to = tt/2 and it is a com- 
bination of Legendre polynomials. For the J = states, 
the dissociation is only possible to continuum levels with 
J = 1, and /(©) is proportional to cos 2 O. For the J = 1 
initial states, the matrix elements of the p —> d-wave 
transitions are larger than the corresponding ones of the 
p —> s-wave transitions. In this case, 1(0) does not have 
zeros but two minima at small but nonzero values which 
are shifted with respect to the two nodes of the second 
order Legendre polynomial. From J = 2 bound states, 
the final continuum level could have J = 1 and J = 3, 
with the transition to J = 1-levels being dominant, and 
1(0) is zero for = tt/2. The quantitatively different 
values of /(©) in both panels illustrate its strong depen- 
dence on the initial bound state due to the overlap of 
the bound state wave function with the highly oscillating 
continuum wave function weighted by the dipole opera- 
tor. For lower lying rovibrational levels, the differential 
cross section is therefore significantly reduced. 
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Fig. 2. Field-free angular distribution for the dissociation process 
to a continuum state with energy E ~ 1.36 mK, which is part of the 
spectrum where the J = 2 shape resonance appears. The initial 
states are (i/, 0, 0) (a) and (z/, 1,0) (b) with vibrational quantum 
numbers v = 40 (full line), 45 (dotted line), 50 (dash-dotted line), 
and 52 (dashed line). 



The effect of increasing the dissociation energy is man- 
ifested in the interplay between the allowed transitions, 
and, thus, in a variation of /(©). In Figs.[2fa) and (b) we 
present the angular distribution from bound states with 
J = and 1, respectively, within the vibrational bands 
v = 40, 45, 50 and 52, for a continuum energy in the 
vicinity of the shape resonance. Note, that a semiloga- 
rithmic scale is used to cover the wide range of variation 
of 1(0). Taking the bound state belonging to the last 
vibrational bands gives rise to variations of /(©) of more 
than two orders of magnitude. If the continuum energy 
is increased /(©) is also increased, as can be seen from 
a comparison of the results for the bound levels (40, 0, 0) 
and (50, 0, 0) presented in Figs.QJa) and (b) and those in 
Fig.[2]^a). Due to the field- free selection rules, the transi- 
tion to the shape resonance is not allowed for the (z/, 0, 0) 
states, whereas it has a dramatic impact on the dissoci- 
ation from the (i/, 1,0) levels, cf. Fig. Efb). Indeed, /(©) 
shows the same dependence on as in Figs. []Ja) and 
(b), but its absolute value is increased by several orders 
of magnitude. For the (50, 1, 0) state, we find 1(0) = 36.8 
and 2.1 x 10 6 for continuum energies 0.3 and 1.36 mK, 
respectively. The resonance wave function has a larger 
probability at short distances, which gives rise to a larger 
overlap with the bound state wave function. In the ab- 
sence of the field, this resonance level is also accessible 
from a bound /-wave, and the corresponding angular dis- 
tribution is also significantly increased compared to other 
spectral regions. 

On exposing the molecule to an additional static elec- 
tric field, the hybridization of the angular motion takes 
place, the selection rule on J does not hold any more, 
and the process of dissociation is altered severely. We fo- 
cus here on the transition from three bound rovibrational 
levels within the field- free vibrational band v = 50, and 
for rotational quantum numbers J = 0, 1, and 2, with 
M = 0, and we consider F < 10.28 KV cm -1 . 

Before discussing the results for the angular distribu- 
tions, it is illuminating to analyze the evolution of the 
above states as the field is varied. The (50, 0, 0) level is a 
high-field-seeker, and for F = 10.28 kV cm -1 , it shows a 
moderate orientation with (cosO) = —0.451. In the con- 
sidered regime of field strengths, the levels (50, 1,0) and 
(50,2,0) are low-field-seekers, and they are not essen- 
tially oriented, (cos 9} = 0.021, and 0.009, respectively, 
for F = 10.28 kV cm -1 . Based on our previous work[32j 
about highly excited rovibrational levels of LiCs in a 
static electric field, we can safely assume that the con- 
sidered field strengths do not yield a significant impact 
on the vibrational dynamics. Hence, in the framework of 
the effective rotor approximation [28|, for a certain field 
strength F, the wave function of the (y, J, M) state can 
be written as 

N-l 

^ vJM (R,0^) « ^o,o(fl) Cf M Y jM (0^), (5) 

j=o 

where Yjm(0, (j>) is the normalized-to-unity spherical har- 
monics with indices j and M, and ifj v o y o(R) are the 



field- free vibrational wave function of the state (v,0,0). 
The Cj M provide the contribution of the partial wave 
(j, M) to the rotational dynamics. They depend on the 
field strength, and satisfy the normalization condition 
J2j I Cj M | 2 = 1 . Note that the above sum should contain 
an infinite number of terms, but for reasons of address- 
ability we assume that only the first N partial waves are 
needed to properly describe the corresponding dynam- 
ics. The evolution of the weights |C^| 2 , j = 0, . . . , 3 for 
the states (50, J, 0) with J = 0, 1 and 2 are presented for 
varying F in Figs.[3fa), (b), and (c), respectively. For the 
(50, J, 0) level, |Cj | 2 remains roughly constant and close 
to 1 for F < 10 kV cm -1 , decreasing thereafter. More- 
over, its contribution is dominant for F < 51.4, 25.7 and 
30.8 kV cm -1 and J = 0, 1 and 2, respectively. The 
weights of the other partial waves, |C^| 2 with j ^ J, 
monotonically increase with increasing F. 
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Fig. 4. Field-dressed angular distribution for the dissociation pro- 
cess to a continuum state with energy E ~ 0.3 mK (a) and 1.36 mK 
(b), the initial state being (50, 0, 0), and F = (long-dashed line), 
0.514 (dash-dotted line), 2.57 (short-dashed line), 5.14 (full line), 
and 10.28 kV cm -1 (dotted line). 
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Fig. 3. Contribution of the partial s (*), p (□), d (•), and / (x) 
waves to the wavefunctions of the states (50, 0, 0) (a), (50, 1, 0) (b), 
and (50, 2, 0) (c) as a function of the strength of the static electric 
field. 

Figures |4ja, b), [5fa, b), and [6fa, b) show I(Q) for 
the two continuum energies 0.3 and 1.36 mK, and for the 
initial levels (50, J, 0) with J = 0, 1, and 2 as well as field 
strengths F = 0, 0.514, 2.57, 5.14 and 10.28 kV cm" 1 . 
For the energy 0.3 mK, the dissociation process is weakly 
affected by the field, see Figs, ffl(a), [5fa) and[5](a). For 
the initial bound state (50,0,0), I(Q) looses its cos 2 6- 
shape as the field is increased, and it shows significantly 
larger values for > 7r/2, compared to < tt/2 for 
F > 2.57 kV cm -1 , which is due to antiorientation. In 
particular, we have 1(0) = 35.6 and 42.4 and I(tt) = 35.6 
and 49.0, for F = and F = 10.28 kV cm -1 , respectively. 
For the (50,1,0) and (50,2,0) states, this effect is only 
appreciable for F = 10.28 kV cm -1 . 

If the dissociation takes place to the spectral region 



close to the shape resonance, the additional static field 
provokes a giant enhancement of the angular distribu- 
tion of the J = and 2 bound states, see Figs. Hl(b), 
and [51(b). For the initial state (50,0,0), 1(0) increases 
by three orders of magnitude when the field is increased 
from F = to 10.28 kV cm -1 , and the corresponding 
value for the angular distribution at zero are 1(0) = 250 
to 3.2 x 10 5 , respectively. The impact on the dissoci- 
ation behaviour of the (50,2,0) state is some what re- 
duced: we obtain 1(0) = 207 and 1.6 x 10 4 for F = and 
10.28 kV cm -1 , respectively. The large enhancements of 
the cross sections and angular distributions are due to 
the field-induced admixtures of p- wave (J = 1) contribu- 
tions to the (50, 0, 0) and (50, 2, 0) states which couple 
directly to the d-wave (J = 2) shape resonance. Never- 
theless, for F = 10.28 kV cm -1 , the weight of the p-wave 
to the (50, 0, 0) and (50, 2, 0) wave functions are still rela- 
tively small, with |C^| 2 = 0.164 and 0.051, respectively. 
Note that these large enhancement of the angular dis- 
tribution are obtained for field strengths which are well 
within experimental reach. Furthermore, for each vibra- 
tional band, exits a field strength that gives rise to the 
proper hybridization of the angular motion and conse- 
quently to an enormous increase of the the dissociation 
probability for the rotational states J / 1 or 3. The 
angular distribution looses its reflection symmetry with 
respect to tt/2 due to the large contribution of the tran- 
sition to the shape resonance. Analogously to the case of 
dissociation to the continuum with energy E = 0.3 mK, 
the angular distribution emerging from the (50, 1, 0) state 
is very weakly affected by the external field. The field- 
dressed /(O) retains its large values, and only a minor 
asymmetry is observed for strong fields. 



IV. CONCLUSIONS 

The dissociation of polar molecules in their electronic 
ground state by the absorption of a single photon to 
produce two ground state atoms has been investigated 
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Fig. 5. Same as in Fig. |4]but for the initial state (50, 1, 0). 
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Fig. 6. Same as in Fig. [4] but for the initial state (50,2,0). 



in the presence of an additional electrostatic field, tak- 
ing as a prototype the LiCs dimer. We have performed 
a fully rovibrational description of the influence of the 
electrostatic field on the molecule, whereas the interac- 
tion with the laser field has been treated by perturbation 
theory. We concentrate on the cold regime with contin- 
uum energies around 1 mK, and on bound levels being 
vibrationally highly excited, 40 < v < 52, and rotation- 
ally cold, J = 0,1,2 and M = 0. Our study focused 
in the experimentally accessible range of field strengths 
F = 0.514-10.28 kV cm" 1 . 

In the absence of the electrostatic field, the dissociation 



on the cold regime which includes several continuum par- 
tial waves such as p — >• s or p —> d transitions. Hence, the 
angular distribution of the states (V, J, 0), with J = 0, 1, 
and 2 shows values of the same order of magnitude. The 
field- free 1(0) is a linear combination of at most two 
Legendre polynomials, being reflection symmetric with 
6. A strong dependence on both the bound and contin- 
uum states is encountered. The presence of a shape res- 
onance gives rise to an enhancement of several orders of 
magnitude for the corresponding cross section of those al- 
lowed transitions. By turning on an electrostatic field the 
dissociation process is drastically altered due to the hy- 
bridization of the angular motion. The transitions from 
the states emerging from the field-free J = and 2 states 
to the shape resonance become allowed and for moderate 
fields, the corresponding angular distribution is several 
orders of magnitude larger than its field- free counterpart. 
Specially, field strengths of a few kV cm -1 provoke gi- 
ant enhancements of 1(0) for the initial states (z^, 0,0), 
whereas a few tens of kV cm -1 are needed to observe a 
similar impact on 1(0) if the process starts from (^, 2, 0) 
levels. 

Although our study focuses on the LiCs dimer and on 
the direct dissociation taking place within its electronic 
ground state, we stress that the above-observed physical 
phenomena are expected to occur equally for other po- 
lar molecules. Shape resonances are typical features in 
the continuum, and their rotational quantum numbers 
depend on the bound states of the highly excited vibra- 
tional bands giving plenty of freedom to the occurrence 
of the observed ehancement effect. 
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